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ABSTRACT 

Since the discovery of the first short-hard 7-ray burst afterglows in 2005, the handful of observed 
events have been found to be embedded in nearby (z < 1), bright underlying galaxies. We present 
multiwavelength observations of the short-duration burst GRB 060121, which is the first observed 
to clearly outshine its host galaxy (by a factor > 10 2 ). A photometric redshift for this event places 
the progenitor at a most probable redshift of z = 4.6, with a less probable scenario of z = 1.7. In 
either case, GRB 060121 could be the farthermost short-duration GRB detected to date and implies 
an isotropic-equivalent energy release in gamma-rays comparable to that seen in long-duration bursts. 
We discuss the implications of the released energy on the nature of the progenitor. These results 
suggest that GRB 060121 may belong to a family of energetic short-duration events, lying at z > 
1 and whose optical afterglows would outshine their host galaxies, unlike the first short-duration 
GRBs observed in 2005. The possibility of GRB 060121 being an intermediate duration burst is also 
discussed. 

Subject headings: gamma rays: bursts; gamma-ray bursts: individual (GRB 060121) 



1. INTRODUCTION 

Since 1993 7-ray bursts (GRBs) GRBs have been clas- 
sified into two subgroups according to the observed du- 
ration and hardness-ratio derived from their gamma- 
ray spe ctra: short-hard (SGRB ) and long-soft (LGRB) 
events ijKouveliotou et al] ll993). 

The first X-r ay det ections of SGRB afterglows by Swift 
llOehrels et alJl200l and HETE-2 (IHiorth et alJ 12003 
lvnfasen^cT^t*aiTl2005|) in combination with optical and 
radio detections have su ggested that the y release less en- 
ergy than the LGRBs i|Fox et al J 12005]) (E 7 (SGRB) ~ 
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10 49 erg vs. E 7 (LGRB) - 10 51 erg) and that they may 
originate from the coalescence of neutron star (NS) or 
NS-black hole (BH) binaries at cosmological distances. 
Until GRB 060121, the detected optical afterglows of 
SGRBs have been comparable in brightness to their host 
galaxies or dimmer, unlike the case of long-soft bursts 
that may outshine their hos ts by a factor of up to ~ 
10 6 (jvan Paradiis et aT1l2000lL The distance scale is also 
considered to be different, having all short bursts with 
definite rcdshifts been zsgrb < 0.6^ while long bursts 
show a broader range, mainly between 0.5 and 5, with 
< zlgrb >^ 1.8. 

2. OBSERVATIONS AND DATA REDUCTION 

GRB 060121 was detected at 22:24:54.50 UT on 21 
January 2006 by FREGATE, W XM and SXC instru - 
ments aboard the HETE-2 mission l)Arimoto et al.l2006[h 
The position of the gamma-ray event was distributed us- 
ing the GRB Coordinates Network (GCN) 13 s later 
The s pectral peak energy Ep = 120 ± 7 keV (Boe r et alJ 
2006), together with a duration of 1.97 ± 0.06 s initially 
classified it in the SGRB group of events (Fig. 1). It was 
also observed by Konus/ WIND llGolenets kii et alf2006a) 
and followed up by Swift /XRT (|Mangano et all 12006 L 
whose observations substantially improved the initial 28' 
HETE-2 error bo x radius and helped t o identify the op- 
tical counterpart l)Malesani et al.ll2006D . 

On receipt of the initial alert by HETE-2, a mosaic 
of images was triggered at the 1.5 m telescope of Ob- 
servatorio de Sierra Nevada (OSN) in order to map the 
entire error box. The first detection of the afterglow 
was obtained only 22 minutes after the 7-ray event. 
Complementary observations were requested with the 2.2 

t We point out that there is still an ongoing debate about the 
redshift of G RB 050813, which could lie at a redshift of z ~ 1.8 
IBergerl2006l) . although no optical or radio counterpart was found 
to support this scenario. 
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Fig. 1. — Burst duration vs. spectral hardness diagram. Large 
black circles are the locations of the three currently established 
HETE-2 short GRBs superimposed on the distribution of 1973 
BATSE short and long GRBs (small grey dots). 







TABLE 1 






Observations of the afterglow of GRB 060121 


Mean Date 


Band 


Tel. 


Int. time 


Flux 


(Jan 2006) 






(s) 


0* Jy) 


22.2443 


K 


4.2m WHT 


750 


17.1±1.4 


23.2402 


K 


4.2mWHT 


1000 


6.3±1.6 


23.3308 


K 


3.5mARC 


3600 


7.48±0.65 


27.2588 


K 


3.5mARC 


3600 


<2.13 


21.9495 


I 


1.5mOSN 


120 


19.3±4.4 


21.9633 


1 


1.5mOSN 


120 


10.0±2.8 


22.0458 


1 


1.5mOSN 


300 


10.8±2.5 


22.1162 


1 


1.5mOSN 


6x300 


4.55±0.83 


22.2539 


I 


1.5mOSN 


5x300 


<2.5 


22.0493 


R 


2.2mCAHA 


600 


3.70±0.97 


22.0963 


R 


2.2mCAHA 


2x600 


1.14±0.33 


22.1590 


R 


2.2mCAHA 


2x600 


1.23±0.49 


22.2385 


R 


2.2mCAHA 


3x600 


1.35±0.37 


23.1828 


R 


1.5mOSN 


12x900 


<0.8 


24.0804 


R 


2.2mCAHA 


6x900 


0.55±0.14 


22.0885 


V 


2.2mCAHA 


2x600 


<1.1 


22.1952 


V 


2.2mCAHA 


5X600 


<0.7 


22.0336 


B 


2.2mCAHA 


600 


<1.7 


22.1560 


B 


2.2mCAHA 


7x600 


<0.9 


22.0258 


U 


2.2mCAHA 


600 


<3.3 


22.1200 


U 


2.2mCAHA 


5x600 


<2.1 



m telescope of the German-Spanish Calar Alto Obser- 
vatory (CAHA), the 4.2 m William Herschel telescope 
(WHT) at Roque de los Muchachos Observatory and the 
3.5 m Astrophysical Research Consortium (ARC) tele- 
scope at Apache Point Observatory. Equatorial coordi- 
nates of the optical near-infrared (nIR) afterglow yielded: 

R.A.(J2000)=09 ,l 09 m 52?02, Dec.(J2000)=+45°39'45."9 
(0. 5 uncertainty at la level). 

Optical images have been photometrically 

calibrated against Sloa n Digital Sky Survey 

ijAdelman-McCarthv et al.1 120061) applyin g the cor- 
responding transformations ijJester et al.1 l2005|) for 
our photometric system. For nIR images we have 
used field stars from the 2 M icron All Sky Survey 
catalogue llSkrutskie et aLl l2~006). A Galactic extinction 
l|Schlegel et alJ 11998(1 correction of E(B-V)=0.014 is 
applied and magnitudes hav e been converted to flux 
densi ty units (Jy) for clarity ifFukueita et al.llT995t iCoxl 
2000). The photometric data are displayed in Table 1. 

The X-ray light curve and spectrum were obtained 
from the Swift/XRT data. The spectrum was fitted 
with a power-law plus a fixed Galactic hydrogen col- 
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Fig. 2. — Spectral flux distribution. A fit of the data points 
with an extincted power law and a Lyman-a break returns two 
peaks of probability. The main one, covering a probability of 63% 
is centred at z = 4.6 with rest frame extinction of Ay = 0.3 ± 0.2 
(top) and the secondary with an 35% is centred at z = 1.7 with 
A v = 1.1 ± 0.2 (bottom). 

umii density (N^q^ = 1.7 x 10 20 cm -2 ) and an in- 
trinsic column density iV H ( int ) at a varying redshift in 
the range 0.1 < z < 6.0 . A spectral index (F ~ v~ & ) 
of = l.lOtg;^ (x 2 /d.o.f. = 36/29) was derived from 
the fit with an intrinsic column density that, depending 
on the selected redshift scenario, ranges from Ah( z= i.5) 

cm" 2 to AT H ( Z= 4.6) 



0.39™ xlO 22 



2 to AW. m = 2.9±;' 3 xl0 22 c 



3. RESULTS 



The detections in I, R and K bands and the upper 
limits imposed for the U, B and V bands, allowed us to 
construct a spectral flux distribution (SFD) at a mean 
epoch of 2.5 hours after the burst (Fig. 2). The nIR K 
band point is extrapolated from a near epoch using as ref- 
erence a quasi-simultaneous R passband detection. The 
data were fit with a power l aw spectrum, a superposed 
intrinsi c extinction llPeil992D and a Lyman-a blanketing 
model (|Madaull 19951) at a varying redshift. The slope of 
the powerlaw was chosen to be (3 opt = 0.60 ± 0.09, as de- 
rived from the X-ray spectra, assuming i/ opt < v c < vx at 
a pre-break e poch, following the standard fireball model 
prescription (jSari et alJl!999|) . This is confirmed by the 
modelling described bellow. We obtained two probabil- 
ity peaks in our redshift study. The main one (with a 
63% likelihood) places the burst at z = 4.6 ± 0.5 with an 
intrinsic extinction of Ay = 0.3 ±0.2 magnitudes. A sec- 
ondary peak (with a 35% likelihood) would imply that 
the afterglow lies at a z = 1.7 ± 0.4 and A v = 1.4 ± 0.4 
magnitudes. In either case, GRB 060121 is the farther- 
most short duration GRB detected to date. A redshift of 
z < 0.5 has a likelihood < 0.5% and implies extinctions 
Ay > 2.0. Adopting a flat cosmology with = 0.73, 
Qm = 0.27 and Hq — 71 km s" 1 Mpc^ 1 and consid- 

we 



ering a 7-ray fluence of 4.71^q'^ x 10 



° erg cm 



obtain an isotropic- equivalent energy release in 7-rays 
of E jAso = 2.18™ x 10 53 erg (£ 7 , lso = 4.30±° j° x 10 52 
erg) for a redshift of z = 4.6 (1.7), 2(1) orders of magni- 
tude higher than other E 7t i so values determined for pre- 
vious SGRBs (Table 2 ). This is compar able to the values 
measured for LGRBs (|Frail et al.ll20QjJ) . 
We have modelled the GRB 060121 afterglow follow- 
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TABLE 2 

Physical properties of the short 7-ray bursts detected to date 



GRB Rcdshift z T90 (s) Flucncc (erg cm 2 ) E 7j j so (erg) Lx (erg s 1 ) Fa/g' 



050509B 


0.225 


0.04 


9.5 x 10" 9 


4.5 x 10 48 


< 7 x 10 41 


< 0.005 


050709 


0.160 


0.07 


2.9 x 10~ 7 


6.9 x 10 49 


3 x 10 42 


~ 1.0 


050724 


0.258 


3.00 


6.3 x 10~ 7 


4.0 x 10 50 


8 x 10 43 


~ 0.2 


050813+ 


0.722/1.8? 


0.60 


1.2 x 10~ 7 


6.5 x 10 50 


9 x 10 43 


< 0.15 


051221A 


0.546 


1.40 


3.2 x 10" 6 


2.4 x 10 51 


6 x 10 44 


~ 1.0 


060121 


1.5 


1.97 


4.7 x 10~ 6 


2.9 x 10 52 


8 x 10 45 


~ 20.0 




4.6 






2.1 x 10 53 


6 x 10 46 




060313 


< 1-7 


0.70 


1.4 x 10~ 5 


< 1 x 10 53 




> 3.0 



Note. — The table displays, in columns: Name of the burst, redshift, duration of the gamma ray emission, measured gamma-ray fluence, 
isotropic- equivalent 7-ray energy, isotropic-equivalent luminosity observed in X-rays 10 hours after the burst and the fr action of afterglow 
flux 12 hours after the b urst and the host galaxy flux (both in R band). The compilation is based on this work and IFox et al.l (2005); 
ISchadv fc Paganil <2006l) ; iGolcnctski TeTan 12006b): Sod erberg eTafl 12003) and references therein. GRB 050813 values are calculated for 
z=0.722. 
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Fig. 3. — Light curve of the afterglow in the near infrared (K), 
visible (R & I) and X-ray bands. The figure shows the result of 
the fit of the model in the most probable high-z (4.6) case which 
has been optimized with values of p = 2.06, 9o = 0.°6 and n = 0.1 
cm -3 and in the low-z (1.5) case, which gives a slightly worse fit 
with p = 2.05, 6» = 2.°3 and n = 0.04 cm -3 . Filled symbols are 
data presented in this article whereas empty ones are data from 
the literature. 

ing the prescription of Uohannesson et al.1 (|200(tT ) on the 
basis of the standard fireball model for which 2 energy- 
injections have been included at 0.035 and 0.23 days af- 
ter the onset of the burst. These energy injections are 
required in order to explain the bumpy behaviour seen 
during the first hours (Fig. 3); s imilar features have 
been seen in long GRB light curves ijCastro-Tirado et al.l 
120061) at z ~ 4. The model has been fitted for both 
z = 4.6 and z = 1.7 scenarios with slightly better re- 
sults for the high-z case. Independent on the redshift, 
the multiwavelength model points to a narrow jet with 
half opening angle 8 < 10° in a low density environment 
(10~ 3 cm~ 3 < n < 0.1 cm~ 3 ; the lower density limit fits 
the observations better for the case where the GRB took 
place at z — 4.6 while the upper density value accom- 
modates better the data at z = 1.7), with efficiencies of 
conversion of kinetic to gamma-ray energy 77^ < 0.05. 

4. DISCUSSION 

GRB 060121 has a T 90 duration of 1.97 ± 0.06 s in 
the 85-400 keV energy band and a hardness ratio (HR) 
S B (100-300 keV)/S£; (50-100 keV) = 1.48 ± 0.18. We 



note that the intrinsic duration of GRB 060121 would 
be 1.6 s with HR of 1.7 if placed at z ~ 0.2 (similar to 
GRB 050709), 0.7 s and 3.3 (if z = 1.7) or 0.4 s and 
4.6 (if z — 4.6). We have studied the classification of the 
progenitor of this GRB (as "old" or "young") considering 
eight criteria: (a) duration, (b) pulse widths, (c) spectral 
hardness, (d) spectral lag, (e) energy radiated in 7-rays, 
(f ) existence of a long, soft bump following the burst, (g) 
location of the burst in the galaxy and (h) the type of 
host galaxy. Four (a, e, f, g) criteria provide strong evi- 
dence that GRB 060121 had an "old" progenitor, while 
four (b, c, d, h) are inconclusive. Thus, we can make a 
strong claim, although not conclusive, that GRB 060121 
had an " old " progenitor. Further details on this analysis 
are given in lDonaghv et all l|2006f) . 

If we c onsider the progenit or to be a me rger o f compact 
objects ijEichler et all 119891: lAlov et alJl2005() . the re- 
leased energy needs either a large conversion efficiency of 
the accreted mass into neutrino emission ( ^ 0.05), a large 
accretion disk mass ( J> .1 Mm) or an appropriat e com- 
bination of both factors ijOechslin fc Jank a 2006) . Con- 
sidering a merger at a redshift z > 1.5, there is a higher 
consistency with the theoretical models where the rate 
of NS+NS or NS+BH merg ers follows the star formation 
rate with delays of ~ 1 Gyr ijJanka et all2005|) than with 
those of very old popula tions of NS+NS progenitor sys- 
tems ijNakar et al.l f2005). Furthermore, the high extinc- 
tion derived from the SFD fit indicates that the merger 
probably took place within the galaxy rather than in the 
outer halo or intergalactic medium and thus, the system 
received a small natal kick, although the density of the 
local event environment is relatively low (n ~ 0.1 cm~ 3 ) 
as indicated by the afterglow fits discussed here. Alterna- 
tively, if the energy was extracted vi a a Blandford-Znajek 
process llBlandford &: ZnaiekHl977j) . either the value of 
the dimensionless angular momentum of the central BH 
is a > 0.3, the magnetic field surrounding the BH is 
B <; 10 16 G or the BH has a mass larger than 3M (or a 
combination of these parameters). 

Late observations by the Hubble Space Telescope (HST) 
have shown no trace of the afterglow down to magnitude 
R ~ 28 about 37 days after the bur st but do show an 
underlying galaxy ijLevan et al.l l2006). We estimate the 
probability of finding such a galaxy within the 0. 5 ra- 
dius erro r box (deriv e d from our astrometry) as 1.3x 10 -2 
following Piro et all l)2002t) . A chance association is un- 
likely but not completely negligible. We have reanalyzed 
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Fig. 4. — Probability distribution vs. redshift. The top panncl 
shows the probability distribution for the GRB afterglow, the host 
galaxy and a combination of both which favours the high redshift 
scenario. In the bottom pannel we have used an extinction-based 
prior which cancels the low redshif peak. 

the p hotometry of the galaxy using ColorPro ijCoe et al.l 
2006) obtaining F606W^ A b = 27.35 ± 0.16 (HST wide 
band filter centred at 606 nm) and F160Wab = 24.05 ± 
0.38 (HST wide band filter centred at 1600 nm). Al- 
though the information is very limited we have used this 
photometry to study the probability distribution for the 
redshift of the galaxy us i ng a bayesian photomeric re dhift 
as described bv lBenitezl l)2000D : iBemtez et alJ l|2004ft . We 
obtain two peaks of probabiliy at z = 1.0 and z = 5.2, the 
latter being more probable. By multiplying this proba- 
bility by the one obtained for the afterglow we can assign 
a probability of 70% for the higher redshift case and a 



28% for the lower redshift scenario (see Fig. 4). Applying 
a prior based on w ell determined extinctions for LGRBs 
ijKann et a l. 2006) the likelihood of a high redshift event 
would rise to 98%. This is a soft constrain as it favours 
low extinctions, like it could be expected from the low 
density derived from our model. 

5. CONCLUSION 

These results suggest that assuming that GRB 060121 
were a SGRB, there exists an emerging population of 
short events located at high redshifts and with energies 
comparable to those of long e vents. In this group we 
may also include G RB 060313 (iSchady fc Paganill2006j) 
or GRB 00030 1C (|Jensen et al.l 12001^ . This popula- 
tion would produce afterglows which significantly out- 
shine their host galaxies, with isotropic energy releases 
of -E 7! j so ~ 10 52-53 erg similar to the values observed in 
long events. Furtherm ore, following the classification by 
iHorvath et all 1)20061) G RB 060121 could be classified in 
the in termediate group l)Horvathlll99 8: Balas tegui et alJ 
12001ft of events with a 68% probability (28% for short- 
burst and 4% for long-burst). The relationship between 
this second population of short bursts and the intermedi- 
ate population of GRBs will be determined or excluded 
by future observations. 

We thank the generous allocation of observing time 
by different Time Allocation Committees. We acknowl- 
edge the Swift team for making public XRT data and S. 
Barthelmy for maintaining the GRB Circular Network 
that distributes the GRB alerts. We also acknowledge 
fruitful dicussions with M. Cervino. 
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